Robust control pulses design for electron shuttling in solid state devices 
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In this paper we study robust pulse design for electron shuttling in solid state devices. This is 
crucial for many practical applications of coherent quantum mechanical systems. Our objective is 
to design control pulses that can transport an electron along a chain of donors, and also make this 
process robust to parameter uncertainties. We formulate it as a set of optimal control problems on 
the special unitary group SU(n), and derive explicit expressions for the gradients of the aggregate 
transfer fidelity. Numerical results for a donor chain of ionized phosphorus atoms in bulk silicon 
demonstrate the efficacy of our algorithm. 



I. INTRODUCTION 

Recent years have witnessed the rapid advance of solid 
state devices that take full advantage of coherent quan- 
tum mechanical properties 1 5]. One particular appli- 
cation of such devices is quantum computation, which 
has attracted intensive research interest over the past 15 
years. To implement these devices in practical applica- 
tions, a central task is to generate a quantum state trans- 
fer. For example, it is often desired to transfer the popu- 
lation between different energy levels or different spatial 
locations so that quantum information can be circulated 
and processed on a large scale. 

In this paper we consider the design of robust control 
pulses for coherent electron shuttling in solid state de- 
vices with a one-dimensional array of quantum dots or 
donors. There are a number of proposals which use the 
donated electron of Group V dopants such as phosphorus 
in silicon as qubits [lJ-H- Such dopants can be inserted 
using ion implantation [fj-|8[. Recently developed tech- 
niques using a scanning tunneling microscope have fur- 
ther allowed them to be placed with high precision Q . 
In order to couple donors and achieve greater scalability 
of qubit array size, it is desirable to be able to move 
qubits robustly between distant physical locations 0. 
The objective is to transport the electron along a chain 
of donors so that the encoded quantum information can 
be communicated between distant qubits. Specifically, 
at the beginning of this procedure, an electron is local- 
ized at one end of the chain. Then by applying some 
appropriate external control fields, we seek to shuttle the 
electron to the other end of the chain. Depending on the 
specific physical implementation, the control fields can 
be gate voltage [l(3, EH or tunable on-site energy |l2j . 

To achieve such electron shuttling, Greentree et al. [l(| 
have proposed to use a solid state version of the well- 
known Stimulated Raman Adiabatic Passage (STIR AP) 
for population transfer in quantum optics [13l Il4| . In this 
adaptation to the solid state, which is known as Coher- 
ent Tunneling Adiabatic Passage (CTAP), two Gaussian 
pulses arc applied in a counter-intuitive sequence to real- 



ize the population transfer in an adiabatic manner, i.e., 
starting in an eigenstatc of the system Hamiltonian and 
changing the Hamiltonian sufficiently slowly so that the 
system will remain in the corresponding eigenstatc during 
the entire transfer process. The amplitudes, peak times, 
and standard deviations of the two CTAP pulses have to 
be carefully tuned. We have recently shown that adia- 
baticity is, however, not a requisite condition for achiev- 
ing high fidelity electron shuttling between spatial loca- 
tions. In Ref. [r|, we have applied Lie-Poisson reduc- 
tion to develop a geometric control approach to remove 
the adiabatic condition and to accomplish the quantum 
state transfer with complete fidelity. 

In the current paper we are interested in the design of 
control pulses for electron shuttling that are robust with 
respect to relevant experimental parameters. In many ex- 
periments it is inevitable that some physical parameters 
are not precisely known although we may have confidence 
that they lie in a certain range. This makes it particularly 
important to design control pulses in a robust manner so 
that the electron shuttling process is insensitive to these 
parameter uncertainties. We formulate this here as an 
optimal control problem on the special unitary Lie group. 
We then discretize the uncertainty range and obtain a fi- 
nite collection of state transfer problems, each of which 
takes a different value of the uncertainty parameter. The 
gradients of the aggregate fidelity with respect to these 
control fields are then derived in an analytic form, which 
allows for efficient implementations of gradient types of 
optimization algorithms. We demonstrate the efficiency 
of our algorithm here by numerical studies with realis- 
tic physical parameters relevant to the electron shuttling 
between phosphorus dopant ions in silicon. 



II. PROBLEM FORMULATION 

In this section we provide a general mathematical de- 
scription for electron shuttling in solid state devices, to- 
gether with the key associated mathematical background. 

The underlying physics and potential applications of 
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solid state devices with qubits have been widely discussed 
in the physical community. See, e.g., Refs. [10l - fl2| . For 
a complete quantum description of the system under re- 
alistic conditions, it is necessary to employ the density 
operator p, which is a Hermitian matrix with unit trace. 
The diagonal elements of the density operator correspond 
to the electron populations on each site. The dynamics of 
the density operator is determined by the Liouville-von 
Neumann equation: 



and the right end of the chain as 



P = ~[iH,p], 



(1) 



where J? is a traceless Hermitian matrix which is termed 
the system Hamiltonian. To be specific, we focus here 
on a triple donor system, but note that the development 
and solution shown here can be easily extended to devices 
with more donors. In this case, the term iH is defined on 
the Lie algebra su[3), i.e., all the 3x3 skew-Hermitian 
matrices. In Ref. [lOj, an electron is moved between ends 
of a chain of ionized phosphorus dopants, for which the 
Hamiltonian is given by (setting h — 1): 
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Here A is the energy difference between eigenstates, and 
f2i2 and f^23 arc the coherent tunneling amplitudes be- 
tween eigenstates. 

Define a basis for su(3) as 
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With a rearrangement of order, this choice of su(3) basis 
is seen to be equivalent to the Gell-Mann matrices fl6j |. 
In this basis, the Hamiltonian in Eq. ([5]) can be repre- 
sented as 



iH = -O12X1 
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(4) 



where I3 is the 3x3 identity matrix. We can drop the 
term y/3 since it commutes with all the other terms and 
thus contributes only a global phase. 

Without loss of generality, let us denote the spatial 
state of the left end of the chain as 



Pi 
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The electron shuttling can now be formulated as a steer- 
ing problem, that is, for the dynamical system of Eq. (JTJ), 
we will apply coherent tunneling amplitudes f2i2 and ^23 
as control fields to transfer the density matrix p from the 
initial state pi at the initial time t — to the final state 
Pt at the terminal time t = T. 

For a fixed energy difference A, this problem has been 
solved by the same authors in 15]. In that work we 
developed an efficient numerical algorithm by using the 
Lie-Poisson reduction theorem. However, as noted above, 
in real experiments, it is often the case that the exact 
value of A cannot be determined precisely, e.g., due to 
imperfections in engineering implementations. Instead, 
we may only know that the energy difference A lies in a 
range [A* — A c , A* + A e ], where A* is the nominal value 
and A e is the maximum possible error bound. These two 
values are usually available for a specific physical system. 

In the rest of this paper, we will design robust control 
pulses that can achieve the desired spatial state transfer 
regardless of what the true energy difference is in the 
given interval. 



III. ROBUST OPTIMAL CONTROL 
ALGORITHM 

To solve the aforementioned robust state transfer prob- 
lem, we take a number of sampling points in the uncer- 
tainty interval and then form a collection of state transfer 
problems, each of which has a different energy difference. 
We then apply a gradient algorithm to find the optimal 
solution that solves all these problems simultaneously. 

To this end, we take N equally spaced points {A n }^ =1 
in the uncertainty interval [A* — A c , A* + A e ], that is, 

. j. . 2(n-l) A 

A„ = A — A e + — -A e , 

n N — 1 

and Ai = A* - A e , A N = A* + A £ . For each A„, 
we consider a dynamical system with the Liouville-von 
Neumann equation 



Pn 



Pn\ 



where the Hamiltonian H n is given by 

A„ A„ 
iH n = —U12X1 — ^23-^2 2~-^ 7 2^71^ 8 ' 



(7) 



(8) 



We now have a set of N dynamical systems, which are all 
identical except for a different value of A in each case. 

We want to steer all these N dynamical systems from 
the initial condition pj in Eq. ([5]) to the final state px 
in Eq. © . Denote the state trajectory of n-th system as 
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p n . We can formulate the state transfer for this system 
as the following minimization problem: 

minL n = \\p T - p n (T)\\%, (9) 

where the Frobenius norm is defined as 

\\A\\% = Tr AAl (10) 

We then have 

L n = Tr(p T - p n (T))(p T - p„(T)) t 
= Tr prp% + Trp n (T)pl(T) 

- TrpTpl(T) - Tr p n (T)p ] T . 

It is easy to show that p = p^ and Tr p n (T)pl l (T) = 
1, and thus minimizing L n amounts to maximizing the 
following fidelity function 



max J n = Tr p T p n (T). 



(11) 



The robust state transfer can now be formulated as maxi- 
mization of the aggregate fidelity of all the terminal states 
Pn(T): 



N 



N 



max J = y] J n = Tr p T p n (T). 



(12) 



the control fields in Eq. (|13|) take constant values which 
are equal to those on the left boundary t = tk'- 
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From Eq. ([5]), we obtain 
iH n (k) = -^i 2 {k)X l - fl 23 (k)X 2 



(14) 
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(15) 

Since iH n (k) is constant on the interval [tk, tfc+i], we can 
compute its unitary propogator as 



U n (k) 



-iH n (k)At 



(16) 



It follows that the density operator at the final time can 
be calculated as 

Pn(T) = U n {K-l) ■ ■ ■ U n (0) Pj UZ(p) ■ ■ ■ Ul(K-l). (17) 

To realize the desired robust spatial state transfer, we 
now only need to maximize the aggregate fidelity in 
Eq. (Ti"2")) with respect to the expansion coefficients a m , 
b m , c m , and d m in Eq. flT3J). 



A. Discretization of sinusoidal control fields 

As discussed earlier, we use the coherent tunneling am- 
plitudes O12 and f223 as control fields. In real physical ex- 
periments, there usually exist maximum frequency limits 
on the control signals. We therefore express the control 
fields as a finite summation of harmonics: 

M 



£li 2 (t) = oq + [a rn cosmwt + b m sin mut], 



m—l 

M 



(13) 



f^23 (t) = Co + /J [c m cos TTLiot + d rn sin vault] , 



m—1 



where lu = lixjT. Here the expansions are truncated 
at a value M, which can be chosen so that Mui stays 
within the feasible frequency range. In the case when 
M is sufficiently large, Eq. (|13[) can approximate any 
continuous control function. 

For time varying control fields, there is generally no an- 
alytic method to solve the Liouville-von Neumann equa- 
tion Eq. ([7]). To obtain numerical solutions, a common 
practice is to divide the total time duration into a num- 
ber of small time steps and assume that the control func- 
tions are constant within each step. In particular, for a 
given time duration [0, T], divide it into K equal intervals 
{[t k ,t k+1 ]}*-^ of length At = tu+i - t k = T/K, where 
tk = kAt. On each of these intervals [tfc 5 ifc+i] , assume 



B. Gradient derivations 

We want to apply a gradient algorithm to find the max- 
imizing expansion coefficients. To this end, we need to 
calculate the derivatives of the cost function J with re- 
spect to those expansion coefficients. 

For the ease of notation, let 
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(18) 
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Then the control fields in Eq. (|14[) can be rewritten in 
the following vector form: 



nia(fc) 

n 23 {k) 



V, 



q, 



(19) 



where the matrix functions cos(-) and sin(-) are calcu- 
lated element-wise. Define 
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where 1 is a column vector with all entries being 1 . Then 
we have 



fii 2 = Gp, Q23 = Gq. 



(20) 



Now the optimization variables become two vectors p 
and q, both of which lie in R M+1 . We proceed to derive 
the gradients of the aggregate fidelity J with respect to 
p and q. From Eqs. (fT2|) and (|20|) . we have 



dJ 
dp 



dfl\ 2 \ dJ 



N 



dp ) c?Sli2 



dJ n 

dSlv 



Similarly, 
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dJ rp * -v dJ^ 
dq~ = G 



dn 



23 



Next we need to derive 



and 



dn 12 (fe) d0 2 3(fc) 



r. Define 



Pn (k) = U n (k - 1) • • • U n (0)piUi(0) ■ ■ ■ Ul(k - 1), 
A„(fc) = Ul(k) ■ ■ ■ UZ(K - l) PT U n (K - 1) • • • U n (k), 

where U n (k) is defined in Eq. (fl6|) . and fc = 0, . . . , K— 1. 
In addition, define /9 n (0) = and A n (K) = p^. Then 

J„ = Tr A n (iOp«(#) = Tr A„(X - l)p n (K - 1) 
= ■■■ = Tr A„(1K(1) = TrA„(0)p„(0). 

It follows that 

dJ n = dTrA n (k + l)p n (k + l) 
dn 12 (k) dn 12 (k) 
rfTrA n (fc + l)f/„(fc)p n (fc)t/t( fc ) 



dn 12 (k) 



:TrA„(fc + b 
fdU n (k) 



Pn (k)Ul(k) + U n (k) Pn (k) 



dUl(k) 
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Using the following expression for the derivative of a ma- 
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trix exponential [11 



i(H a +vH h )t 
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Substituting Eq. ^ into (gT]), we get 



dJ n 



dCl 12 {k) 
= TrA„(fc + l) 
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= Tr[p„(fc + l),A n (fc + l)] / e- iH ^ T X ie iH ^ T dT. 

Jo 

(23) 

We can further simplify the calculation of Eq. (|2"3"|) . We 
first note that since H n (k) is a Hcrmitian matrix, it can 
be diagonalized as 



H n (k)=T n (k)T n (k)Tl(k), 



(24) 



where 



r il (fc)=diag{ 7l 1 l (fc),7'(fc),7'(fc)} 

f A„ A n + 3g n (k) A„ - 3g„(fc) 
-diagj-^, - , - 



and the unitary matrix T n (k) can be written as 



T„(k) = 



-n 23 (k)/h(k) n 12 (k )/^g n (k)(g n (k) +A n )/ 2 n 12 ( k)/^/g n (k)(g n (k) - A n ) /2 

-V(g n (k) + An)/(2g n (fc)) V(gn( fc) - An)/(2 gn (fc)) 
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and 



Therefore we can write, 



5 „(fc) = VA2+4f2| 3 (fc)+40f 2 (fc), 



Mfc) = ^ni 3 (fc) + n?2(fe)- 
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(25) 
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where denotes the Hadamard product, i.e., element- 
wise product, of two matrices, and the a6-th element of 
^ n (fc) is exp{i(7^(fc) — Jn(k))r}. Now define a matrix 
$ n (fc), whose a6-th element is given by 



f At 



*" 5 (*)dT 

'cxpM(fc)- 7 £(fc))At} 



a*, 



for a b. 
for = 6. 



This allows Eq. ([25]) to be calculated explicitly: 

At 

io (26) 

=T n (*0((T t t(fc)XiT n (fc)) $„(fc))T r t(fc), 

which in turn yields that 



dfl 12 (k) 

= Tr( [p n (fc) , A„ (fc)]T„ (A) ( (Tt (fc) AT^ (A)) *)Tt (fc) ) . 

(27) 

A similar analysis leads to 
dJ n 



dfl 23 (k) 

= Tr([p„(fc),A„(fc)]r„(fc)((r t t(fc)X 2 T„(fc))0$)T J l(fc)). 

(28) 

We have derived closed form formulae for dJ /dp and 
dJ/dq, i.e., the gradients of the aggregate fidelity J with 
respect to the expansion coefficient vectors p and q. It is 
now straightforward to implement gradient types of algo- 
rithms such as the gradient descent algorithm, the con- 
jugate gradient algorithm, or the BFGS algorithm [l8j . 



IV. NUMERICAL OPTIMIZATION 

In this section we apply the gradients derived in the 
preceding section to design the robust control fields that 
can realize the desired population transfer in solid state 
devices. 

We consider the ionized donor chain that was discussed 
in Ref. [uj. Typical values of the energy difference A 
are several meV, while the control fields f2i2 and 
can be varied in the magnitude of 10 -2 meV. Realistic 
parameter values allow us to assume a nominal value for 
A* of 2.72 meV, with the actual value of A deviating from 
the nominal value by up to 20%. We further assume that 
the population transfer needs to be accomplished within 
100 ns, and the maximum feasible frequency for control 
fields is 0.1 GHz. These constraints lead to the control 
field expansions in Eq. (|14j) needing to be truncated at 
M = 10. 

We now discretize the total time duration [0, 100] ns 
into 100 small time steps, each with length 1 ns. Take 
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FIG. 1: Control pulses obtained with no robust design. (A) 
Spatial state transfer when A = A* = 2.72 meV. (B) Control 
pulses: blue solid line corresponds to f2i2, green dashed line 
to f223- (See electronic version for color plots). 



11 evenly distributed sampling points from the uncer- 
tainty range [0.8A*, 1.2A*] meV. Given these parameter 
settings, we can apply a gradient algorithm with fixed 
step size to solve for the optimal control pulses. 

As a reference, we first consider the case with no robust 
design, i.e., optimizing for the point A = A* only. The 
corresponding population transfer and control fields are 
shown in Fig. [TJA) and (B), respectively. To test the ro- 
bustness, we apply these control pulses to all 11 sampling 
points in the uncertainty range [0.8A*, 1.2A*] meV. The 
results of these simulations are plotted in Fig. [5] It is ev- 
ident that when the actual value of A is unknown within 
this range, the electron cannot be successfully transferred 
from left to right, except in the case when (coinciden- 
tally) A = A*. Note that in each of the unsuccessful 
transfers, a full transfer is achieved at some point before 
the final time T. However, the oscillatory nature of the 
populations leads to a reversal of the transfer. There- 
fore, experimentally, a number of different transfer times 
would have to be attempted for a given pulse sequence 
in order to assess the possibility of a complete transfer 
and to determine the optimal time. Furthermore, noise 
in any element of the Hamiltonian may cause the optimal 
transfer time for a given pulse sequence to be different 
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FIG. 2: Robustness test for the control pulses in Fig. [T] spa- 
tial state transfers for 11 evenly distributed A's in the range 
[0.8A*, 1.2A*], where A* = 2.72 meV. 



for each individual experiment. 

Next we apply the robust control pulses design devel- 
oped above. The optimization results for this scheme are 
shown in Fig. [3j The robust controls are seen to be about 
an order of magnitude larger than the controls for A* 
only. Most importantly, it is clear that whatever value 
of energy difference A within the ±20% deviation range 
of the nominal value A* = 2.72 meV is employed, the re- 
sulting robust control fields can transfer the population 
with almost perfect fidelity. The robust controls also have 
the advantage that the populations do not oscillate as in 
Fig. [2j Slight changes in transfer time would therefore 
not affect the population transfer, a useful robustness 
feature from the experimental perspective. Finally, we 
note these pulses also perform well outside the range for 
which they were defined. For example, if the real uncer- 
tainty level is ±25% instead of ±20% in the design, the 
spatial state transfer still has acceptable performance, as 
shown in Fig. |4l 



sign significantly improves the performance of an electron 
shuttling protocol, achieving near perfect state transfer 
across a realistic range of Hamiltonian parameters for a 
phosphorus-doped silicon system. 



0.8 



s 0.6- 



0.4 



0.2 





































Pll 

- "P22 
P33 










/ 











20 



40 60 

Time (ns) 



80 



100 



(A) 




-0.005 



40 60 

Time (ns) 



100 



(13) 



FIG. 3: Design of robust control fields. (A) Spatial trans- 
fer for 11 evenly distributed values of A's in the range 
[0.8A*, 1.2A*], where A* = 2.72 meV. (B) Robust control 
pulses: blue solid line corresponds to S7i2, green dashed line 
to f223- (See electronic version for color plots). 



CONCLUSION 
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